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ABSTRACT 


An  investigation  into  the  general  nature  and  properties  of  phase 
trajectories  of  third  order  linear  feedback  control  systems  has  been 
conducted. 

The  geometry  of  eigenvectors,  eigenplanes,  and  the  conic  surface 
named  the  eigencone  was  investigated  with  respect  to  system  root 
parameters.  An  analysis  has  been  made  of  the  orientation  of  eigenplanes 
of  systems  with  three  real  roots,  and  of  complex  eigenplane  geometry  as 

the  damping  ratio  £,  and  undamped  natural  frequency  uu  were  varied. 

n 

A  method  is  introduced  for  determining  system  roots  that  will  locate  an 
eigenplane  or  complex  eigenplane  in  a  pre-determined  location  in  third 
order  error  space. 

Analog  and  digital  computer  programs  used  to  obtain  phase  trajec¬ 
tories  of  third  order  linear  systems  are  presented.  The  effect  of  root 
location  in  the  s-plane  on  phase  trajectories  is  discussed.  A  thorough 
analysis  of  the  effect  of  initial  conditions  on  real  root  and  complex 
root  phase  trajectories  is  presented  along  with  numerous  photographs  of 
three-dimensional  phase  trajectory  models. 

The  inter-relation  of  phase  trajectories,  eigenvectors,  and  eigen¬ 
planes  is  analyzed.  An  analytical  investigation  of  the  response  of  a 
complex  root  system  to  a  step  input  was  made  as  a  first  step  towards 
possible  application  of  error  space  geometry  to  discontinuous  control 
systems . 
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Chapter  ] 


INTRODUCTION 

In  recent  years  many  papers  have  been  written  with  regard  to  the 
analysis  and  design  of  discontinuous  feedback  control  systems.  In  many 
of  these  papers,  various  schemes  have  been  proposed  for  changing  the 
system  damping  at  some  optimum  time  to  produce  a  deadbeat  response. 

However,  the  rigorous  mathematical  developments  of  the  switching  criteria 
in  some  of  these  proposals  tend  to  obscure  somewhat  the  overall  picture. 
While  this  investigation  does  not  concern  itself  with  discontinuous 
systems,  the  phase  trajectories  of  third  order  linear  feedback  control 
systems  have  a  definite  relationship  to  the  eigenvector  geometry  of 
three-dimensional  error  space  which  lends  itself  to  easy  visualization 
of  switching  surfaces. 

The  general  nature  of  the  problem  was  to  investigate  the  phase 
trajectories  of  third  order  linear  feedback  control  systems  in  response 
to  various  initial  conditions.  It  was  intended  originally  to  devote 
the  major  portion  of  the  investigation  to  phase  trajectories  of  systems 
characterized  by  one  real  root  and  a  pair  of  complex  conjugate  roots 
and  only  to  verify  Han's  work  (Ref.  1)  regarding  the  trajectories  of 
systems  with  three  real  roots.  However,  during  the  course  of  the 
investigation,  it  was  found  that  the  phase  trajectories  of  both  these 
types  of  systems  have  a  definite  relationship  to  the  eigenvectors 
associated  with  the  system  and  to  certain  planes  and  surfaces  defined  by 
these  eigenvectors.  It  was  found  possible  to  describe  any  phase  trajectory 
of  any  third  order  linear  system  in  response  to  any  set  of  initial 
conditions  in  terms  of  their  relationship  to  the  eigenvectors  associated 
with  the  system. 

The  use  of  eigenvectors  to  define  surfaces  and  planes  in  phase 
space  is  not  a  new  concept  but  this  application  has  been  largely  confined 
to  mathematical  definitions  and  manipulations  The  approach  adopted 
during  the  course  of  this  investigation  was  to  construct  various  repre¬ 
sentative  types  of  phase  trajectories  and  eigenvectors  in  a  three-dimen¬ 
sional  coordinate  system  and  observe  them  in  relation  to  each  other. 
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In  the  analysis  and  design  of  second  order  feedback  control  systems, 
in  particular,  those  systems  with  non-linearities,  the  use  of  the  phase 
plane  has  become  wide-spread  It  is  conceivable  that  the  use  of  phase 
space  could  also  become  a  valuable  tool  in  the  analysis  and  design  of 
third  order  and  higher  feedback  control  systems- 

The  basic  concept  of  phase  space  is  the  choosing  of  a  particular  set 

of  time  variables  of  the  control  system  as  the  coordinate  axes  of  this 

space.  In  most  feedback  control  systems,  the  system  error  automatically 

seeks  zero,  and  the  phase  space  coordinates  are  normally  chosen  as  error, 

E,  error  rate  or  velocity,  E,  and  the  (n-l)  derivatives  of  error  necessary 

for  control  system  of  order  n.  The  behavior  of  the  system  is  then 

#  •*  (n-i) 

described  by  the  spatial  coordinates  E,  E,  E,  ...  E.  These  coordinates 
are  the  previously  mentioned  time  variables  of  the  control  system,  and 
then  the  system  may  be  described  mathematically  by, 


For  the  third  order  systems  which  this  paper  is  concerned  with.  Equation 
(1-1)  becomes,  in  more  familiar  form, 


(1-2) 


E  +  aE  +  bE  +  cE  *  0 


Thus  at  any  time,  t,  there  is  a  point  in  space,  called  the  "state" 
point,  which  describes  the  behavior  of  the  system,  and  as  time  varies 
from  t  58  0  to  t  *  °°,  this  state  point  traces  a  path  in  this  three-dimen¬ 
sional  space,  called  the  phase  trajectory  As  time,  t,  approaches  infinity, 
the  system  error  seeks  zero,  and  the  phase  trajectory  approaches  the  origin 
of  this  phase  space. 

Since  this  three-dimensional  space  is  based  on  system  error  and  its 
derivatives,  the  phase  space  will  hereafter  be  referred  to  as  "error 
space" . 

In  previous  investigations  using  phase  space,  for  example,  Bogner's 
(Ref.  2),  standard  matrix  techniques  were  used  to  manipulate  phase  space 
into  principal  coordinate  space.  This  coordinate  system  uses  the  eigen¬ 
vectors  as  coordinate  axes.  This  transformation  is  particularly  useful 
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in  the  analysis  of  higher  order  contactor  servo  systems.  Since  this 
investigation  is  based  primarily  on  physical  interpretation  of  phase 
trajectories,  this  transformation  to  principal  coordinate  axes  was  not 
made . 

By  limiting  the  order  of  the  systems  to  be  investigated  to  three, 
and  by  not  transforming  into  principal  coordinate  space,  physical  repre¬ 
sentation  of  phase  trajectories  in  three-dimensional  error  space  was 
possible.  To  aid  in  this  physical  representation  of  phase  trajectories, 
the  coordinate  planes  and  octants  of  third  order  error  space  are  defined 
as  in  Fig.  1,  and  transparent  lucite  models  of  error  space  similar  to 
Fig.  1  were  constructed 

The  various  phase  trajectories  investigated  were  obtained  from  both 
digital  and  analog  computers.  Wire  models  of  these  trajectories  were 
constructed  in  the  coordinate  system  represented  by  the  plastic  models. 


0 
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Fig.  1  Coordinate  Planes  of  Thi^d  Order  E^ro^  Space 
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Chapter  11 


EIGENVECTORS  IN  THIRD  ORDER  ERROR  SPACE 


General 

It  is  well  known  that  for  each  real  root  of  a  third  order  system 
there  is  an  eigenvector  which  corresponds  to  a  phase  trajectory  that  is 
a  straight  line.  (See  Ref.  3.)  For  a  system  with  three  real  and  distinct 
roots,  r^  <  r^  <  r  there  are  three  eigenvectors.  The  eigenvector 
associated  with  the  smallest  root  r  is  referred  to  as  the  slow  eigenvector 
since  the  time  constant  of  r^  is  larger  than  those  of  the  other  two  roots 
and  produces  a  slower  response.  Similarly,  since  r^  has  the  shortest 
time  constant,  its  associated  eigenvector  is  termed  the  fast  eigenvector. 
Since  the  time  constant  of  r^  is  some  intermediate  value  between  that  of 
r^  and  r^>  the  eigenvector  associated  with  it  is  referred  to  as  the 
intermediate  eigenvector.  In  the  case  of  a  system  defined  by  one  real 
root  and  a  pair  of  complex  conjugate  roots,  due  to  the  presence  of  only 
one  real  root,  there  will  be  only  one  eigenvector. 


Determination  of  Eigenvectors 

The  differential  equation  of  a  linear  third  order  feedback  control 
system  may  be  expressed  in  terms  of  the  system  error,  E,  as 

*E  +  aE  +  bE  +  cE  *  0  (II- 1) 

where,  in  terms  of  the  system  roots  r(,  r^  and  r^,  all  assumed  real 
and  distinct, 


a  *  ri  *  r2  *  r3 

b  *  rlr2  +  rlr3  +  r2r3 


rlr2c3 


This  third  order  differential  equation  may  be  rewritten  as  three 
dependent  first  order  differential  equations  by  first  setting  E  =  E  . 

Then , 

E  =  E  -  E2  (II-2a) 

E  *  E2  -  E3  (ll-2b) 

E  =  E^  =*  -cE1  -bE2  -aE  j  (II-2c) 
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The  set  of  equations  thus  obtained  can  be  written  in  matrix  form  as 

Thus  , 


-  1, 

3. 

- 1 

M  • 

_ 1 

*2 

a 

=  3 

0 

0 

-c 


1 

0 

-b 


L 


(I  I  -  3) 


The  matrix  of  coefficients,  A,  has  an  eigenvalue  X  if  and  only  if 
there  exists  a  vector  E  with  components  ,  Ey  and  E  ^  such  that  AE  -  XE , 
where  multiplication  is  matrix  multiplication  and  E  is  considered  to  be 
a  one  column  matrix.  Such  a  vector  E  is  called  an  eigenvector  of  the 
matrix.  Eigenvectors  exist  satisfying  AE  =  XE  only  when  X  satisfies  the 
determinant  (A  -  Xl)  =  0,  where  I  is  the  identity  matrix. 


det  (A  -  Xl) 


Expanding  the  determinant, 


-X 

0 

-c 


1 

~b 


0 

1 

-  (a+X) 


(II-4) 


X3  4*  aX^  4*  bX  4  c  33  0 


(II-5) 


The  expression  obtained  is  the  characteristic  equation  of  the  matrix, 
not  to  be  confused  with  the  characteristic  equation  of  the  feedback 
control  system.  Any  root  X^  of  this  equation  is  an  eigenvalue,  and 
associated  with  each  eigenvalue  there  is  an  eigenvector.  The  character¬ 
istic  equation  factors  into  the  form 


+  r1)(X2  +  r£)  (X3  4-  r}) 


0 


(11-6) 


Thus , 


-  r . 


-  r . 


To  find  the  eigenvector  associated  with  each  root,  each  root  is  sub¬ 
stituted  in  turn  into  the  matrix  equation  (A  -  X^J)  ^ 


- 1 

1 

H- 

1 

0 

1 

ij*\r 
_ 1 

1 

t 

o  o 

-X. 

1 

-b 

1 

-(a+X1) 

X 

E2 

r3 

—  .J 

(II-7) 
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(II -8a) 


-X.  E 

+  E  , 

=»  0 

l  1 

2 

-A.  e! 

+  e! 

=*  0 

i  2 

3 

— V 

-CE1  ‘ 

bE2  “ 

( a+\ 

(II -8b) 
(I  I -8c) 


To  obtain  a  solution,  any  two  of  the  three  equations  may  be  solved  in 
terms  of  one  of  the  variables  and  an  arbitrary  constant  K  assumed  for 
its  value.  The  easiest  choice  here  is  to  solve  Equations  (II-8a)  and 
(II-8b)  in  terms  of  and  to  assume  E  }  ~  K. 


K_ 

X. 


(II -9a) 


®  X.K 

i 


(II-9b) 


The  values  obtained  for  E^,  E^  and  E.^  may  be  checked  by  substitution  into 
Equation  II-8c. 


The  general  expression  for  an  eigenvector  may  then  be  written  as 


E  =  y—  +  KE2  +  X.  KE^  (II - 10) 

i 

where  E^,  E0  and  E^  may  be  considered  as  unit  vectors  along  their 
respective  axes.  In  terms  of  E,  E,  E  and  r^,  where  is  any  distinct 
real  root  of  the  system,  this  expression  may  be  written  as 

IT  *  - —  E  +  KE  -  r.KE  (11-11) 

-  r  i 

i 


From  this  general  expression  of  an  eigenvector  in  three,  dimensional  error 
space,  it  is  seen  that  an  eigenvector  has  a  fixed  orientation  which  is 
dependent  only  upon  the  location  of  the  associated  real  root.. 


Throughout  the  course  of  the  investigation  when  it  was  desired  to 
determine  the  eigenvectors  of  a  particular  system,  it  was  necessary  to 
locate  only  one  point  on  a  particular  eigenvector  to  orient  it  properly 
in  space.  From  the  general  expression  derived  for  an  eigenvector,  any 
arbitrary  point  on  an  eigenvector  is  given  by  the  coordinates  (E,  E,  E)  31 


~7~ 


1.0  and 


For  example,  assume  a  system  with  real  roots  at  r^  ~  .5,  r^  = 
r^  3  3.0.  Here  the  roots  are  defined  as  positive  when  they  are  located 
in  the  left  half  of  the  complex  s-plane.  By  assuming  a  value  for  K,  the 
coordinates  of  a  point  on  each  of  the  eigenvectors  may  be  determined. 

Let  it  be  assumed  that  K=l.  The  resulting  point  on  each  of  the  eigen¬ 
vectors  are  listed  below. 


Eigenvector 

s  low 

intermediate 

fast 


(A,  A,  E) 

(-2,  1,  -.5) 
(-1,  1,  -1) 
(-0.333,  1,  -3) 


Note  that  all  three  eigenvectors  lie  in  octant  2*.  (Refer  to  Fig.  1)  . 
However,  suppose  that  K  *  -1.  Then 


Eigenvector 

slow 

intermediate 

fast 

These  eigenvectors  all  lie 
shown  in  Fig.  2. 


(E,  E,  E) 

(2,  -1,  .5) 

(1,  -1,  1) 

(0.333,  -1,  3) 

in  octant  4.  A  model  of  this  example  is 


If  the  roots  of  the  system  had  been  chosen  as  negative,  i.e., 
located  in  the  right  half  of  the  complex  s-plane,  then  the  eigenvectors 
would  be  located  in  octants  1  and  3f  .  Thus,  for  a  third  order  system, 
all  possible  eigenvectors  of  the  system  are  restricted  to  octants 
1,  2 1 ,  3'  and  4. 


The  Eigencone 

From  observations  of  eigenvector  location  in  three  dimensional 
models,  it  was  theorized  that  there  is  a  conical  surface  in  three  dimen¬ 
sional  error  space  which  is  the  locus  of  all  possible  eigenvectors  for 
a  third  order  system.  In  order  to  prove  this  theory,  the  coordinates 
of  any  arbitrary  point  (E ,  E,  E)  on  an  eigenvector  are  written  in  para¬ 
metric  form  as 

E  *  E  =  K  E  =  -r.K  (11-12) 

-r.  1 

1 


-8 


9 


(II  - 13a) 


Eliminating  the  parameter  K  and  solving  for  r_.  , 


E  +  r.E  -  0 
1 


r . 
i 


E 

E 


E 


(II- 13b) 


Eliminating  r_.  from  these  equations  leaves 

E2  -  EE  -  0  (11-14) 

which  is  the  equation  of  a  right  elliptical  cone.  Since  direct  sub¬ 
stitution  of  the  parametric  equations  for  any  arbitrary  point  on  an 
eigenvector  satisfies  this  equation  of  a  right  elliptical  cone,  the 
surface  of  the  cone  is  the  locus  of  all  possible  eigenvectors  associated 
with  the  roots  of  a  third  order  system 


This  cone,  which  will  be  referred  to  as  the  eigencone,  can  be 
visualized  as  shown  in  Fig.  3  as  the  surface  formed  by  an  infinite 
number  of  eigenvectors  emanating  radially  from  the  origin  of  the 
coordinate  system.  Figure  3  shows  only  that  portion  of  the  eigencone 
determined  by  the  eigenvectors  of  roots  located  in  the  left  half  of 
the  complex  s-plane.  An  exploded  view  of  the  eigencone  section  located 
in  octant  2*  is  shown  in  Fig.  4  Figure  5  shows  this  same  section  of 
the  eigencone  in  its  proper  location  as  viewed  looking  into  octant  4. 


From  these  photographs  it  is  apparent  that  the  apex  of  the  eigencone 
is  located  at  the  origin  of  error  space  and  is  tangent  to  both  the  E  and 
E  axes.  It  is  symmetrical  about  the  line  E  38  E  in  the  E/E  plane  and  is 
bisected  by  this  plane.  It  will  be  noted  in  Fig  5  that  the  surface  of 
the  eigencone  nearest  the  E  axis  is  labled  SLOW  while  the  surface  nearest 
the  E  axis  is  labled  FAST  The  eigenvectors  in  these  regions  correspond 
to  phase  trajectories  which  have  response  times  characteristic  of  these 
labels.  By  referring  again  to  the  parametric  equations  of  a  point  on  an 

eigenvector  Equa t ion  (ll- 12),  it  will  be  noted  that  for  r^  <  1.0,  the 

eigenvector  will  be  located  in  the  SLOW  region  while  for  >  1.0,  the 

eigenvector  will  be  located  in  the  FAST  region.  When  r^  ~  1.0  the 

eigenvector  is  equidistant  from  both  the  E  and  E  axes.  It  should  be 
pointed  out  here  that  the  classification  of  eigenvectors  into  slow. 
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intermediate  and  fast  1=  not  according  to  the  labels  on  the  eigencone, 

but  is  applied  only  to  eigenvectors  associated  with  the  roots  of  a  given 

system  It  is  incorrect  to  arbitrarily  call  all  eigenvectors  slow  if 

they  are  associated  with  r.  <  1.0  or  fast  if  r.  >  1  0.  For  different 

ii 

systems,  the  slow  eigenvector  of  one  may  have  a  faster  response  than 
the  fast  eigenvector  in  the  other.  The  labeling  of  the  eigencone  into 
slow  and  fast  regions  is  intended  to  indicate  the  relative  speed  of 
response  as  related  to  the  eigenvectors,  and  not  to  establish  any  sharp 
line  of  demarcation. 
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Chapter  111 


EIGENPLANES  IN  THIRD  ORDER  ERROR  SPACE 


General 


It  has  been  shown  (Ref.  1,  4,  5)  that  any  two  eigenvectors  in  error 
space  define  a  plane,  passing  through  the  origin  of  error  space,  called 


an  eigenplane  (or  sometimes  a  "hyperplane") .  The  eigenplane  has  an  equation 


one  order  lower  than  the  order  of  the  characteristic  equation  of  the 
control  system.  The  eigenplane  equation  is  obtained  by  suppressing  the 
root  corresponding  to  the  pertinent  eigenvector  in  the  differential 
equation  of  the  system.  Thus,  in  general, 


(TII-1) 


For  the  third  order  case  with  real  distinct  roots,  r^  <  r^  <  r^. 


36  E  +  (r2+r3^  E  +  r2r3  E  23  0  (r^  suppressed) 

P»2  32  E  +  (r^4r^)  E  +  r^r^  ^  ^  0  (r^  suppressed)  (III-2) 

a  E  +  (r  -Hr,,)  E  +  r.^  E  *  0  (rj  suppressed) 


are  the  defining  equations  for  the  three  eigenplanes,. 

The  eigenplanes  in  third  order  error  space  are  also  defined  by  two 
straight  lines,  the  two  eigenvectors  corresponding  to  the  two  roots 
contained  in  the  eigenplane  equation.  That  is,  the  eigenplane  defined 


by  suppressing  the  smallest  root,  r^ ,  whose  corresponding  eigenvector 


is  the  slow  eigenvector,  contains,  and  can  be  defined  by,  the  intermediate 
and  fast  eigenvectors.  In  a  like  manner,  the  other  two  eigenplanes  are 
defined  by  the  combination  of  slow  and  fast  eigenvectors,  and  by  the 
combination  of  slow  and  intermediate  eigenvectors 

When  referring  to  specific  eigenplanes  hereafter  in  this  paper,  the 
MFirstu  eigenplane  will  mean  that  plane  containing  the  intermediate  and 
fast  eigenvectors  formed  by  suppressing  the  smallest  root,  r^ ;  the 
"Second"  eigenplane  will  mean  that  plane  containing  the  slow  and  fast 
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eigenvectors  formed  by  suppressing  the  root,  r>;  and  the  "Third" 
eigenplane  will  mean  that  plane  containing  the  slow  and  intermediate 
eigenvectors  formed  by  suppressing  the  largest  root,  r ^ 


Eigenplanes  of  Systems  With  Three  Real  Distinct  Roots 

In  a  control  system  having  three  real  distinct  roots,  there  will, 
of  course,  exist  three  eigenplanes  in  error  space.  Since  these  eigen¬ 
planes  are  each  defined  by  two  eigenvectors,  they  should  intersect  the 
eigencone  along  lines  corresponding  to  their  two  defining  eigenvectors. 
This  can  be  proved  analytically  by  a  simultaneous  solution  of  the 
equations  of  the  eigencone  and  of  the  eigenplane  as  shown  below: 

-  2  •• 

Equation  of  eigencone:  E  -  EE  -  0 

Equation  of  eigenplane:  |i  3  E  +  AE  +  BE  *0 

where  A  *  r^  +  r^ 


Solving 


the 


equation 


of  the  eigencone 


for  E  yields: 


This  expression  for  E,  when  substituted  into  the  equation  of  the 


eigenplane,  yields; 


E  +  AE  +  B 


=  0 


or 

E2  +  AEE  +  BE2  *  0 

Now,  assume  E  to  have  any  arbitrary  value,  K.  Hence, 

-2  ••  2 

E  +  AKE  +  BK  =*  0,  and  this  expression,  when  solved  by  the  quadratic 
formula,  yields  for  E, 


E  = 


-/IK  t  V(*K)2-4BK- 
2 
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or  in  terms  of  the  roots,  r ^  and  r^ 


t 


-(r;.r2)K  t'WK^  Ztrr.K'+CK’-^rr,  K: 


-cK-rj.K 

2 


±  UO^~r~K) 


-rt  K ~t\K  ±  hK-r^K) 

~  z 

Thus  :  E  =  -  rz  K  ,  -  r,  K 

P  _  E2  _  K 2  __K_Z  _  pK_ 

E  ~~^K’  -r,  K  r2.  ’  r< 

Hence:  (E,E,E)  -  K ,  “  K  r,  )  and  (-p  >  K,  ~Kr2) 

which  are  the  coordina tes  of  any  point  on  the  two  eigenvectors  defining 
the  eigenplane,  |i. 

If  one  were  to  visualize  all  three  eigenplanes  in  error  space,  a 
volume  would  be  bounded  by  these  three  eigenplanes.  This  volume  would 
be  two  pyramids,  in  octants  4  and  2*  with  their  apexes  at  the  origin 
of  error  space  and  their  edges  corresponding  to  the  three  eigenvectors 
of  the  system.  Fig  6  is  a  pictorial  representation  of  just  the 
pyramid  in  octant  21  defined  by  the  three  eigenvectors  and  bounded  by 
the  three  eigenplanes  of  a  system  with  roots,  r  ®  0  333,  r?  53  10, 
and  r3  =*  3.0 

The  orientation  of  this  pyramid  is  such  that  it  will  always  lie 
within  the  eigencone  and  its  size  and  shape  will  depend  on  the  location 
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of  the  system  roots  on  the  axis  of  reals  If  the  system  roots  are  all 
close  together,  the  volume  of  the  pyramid  will  be  small  as  compared 
to  the  volume  enclosed  by  the  eigencone.  Conversely,  if  r^  is  allowed 
to  approach  zero  and  r^  is  allowed  to  approach  infinity,  with  r^  ^  1.0, 
the  resulting  pyramid  will  almost  completely  fill  the  eigencone. 

If  all  three  system  roots  are  small,  i.e.,  less  than  1.0,  the 
pyramid  will  be  relatively  small  and  oriented  near  the  E  axis;  and, 
conversely,  if  all  three  roots  are  greater  than  10,  the  pyramid  will 
also  be  small,  but  oriented  near  the  E  axis. 

Eigenplanes  of  Systems  Having  One  Real  Root  and  a  Complex  Pair  of  Roots 

For  a  control  system  having  one  real  root  and  two  conjugate  complex 
roots,  it  can  be  readily  seen  from  equations  (TTI-2)  which  define  the 
three  general  eigenplanes  that  only  one  eigenplane  can  exist  for  this 
system.  If  r^  and  r^  are  of  the  form  a  +  jb,  the  eigenplanes  and 
cannot  exist  because  of  the  imaginary  coefficients  of  E  and  E  in 
Equation  (ill-l).  Thus  only  one  eigenplane,  |J>  38  E  +  (r,-frr^)E  +  v ir 3  E  58  0 
and,  as  was  shown  in  the  previous  section,  only  one  eigenvector  with 
coordinates | ,  +  k,  +  k  r^J  will  exist  in  error  space.  This  single 

eigenplane  existing  for  this  type  system  will  subsequently  be 
referred  to  ^  s  the  "complex  eigenplane". 

Closer  inspection  of  the  equation  of  the  complex  eigenplane  given 
above  shows  that  since  the  real  root,  r.j  ,  does  not  appear  in  the  equation 
the  location  of  the  complex  eigenplane  in  error  space  is  completely 
independent  of  the  magnitude  of  the  real  root  of  the  system. 

An  investigation  of  the  effects  of  varying  system  parameters  and 
complex  root  location  was  conducted  to  gain  more  knowledge  about  the 
orientation  of  the  complex  eigenplane  in  error  space  In  this  phase  of 
the  investigation,  the  three-dimensional  model  was  invaluable.  By 
setting  any  one  of  the  three  coordinates  of  Equations  (111-2)  equal 
to  zero,  the  traces  in  the  three  coordinate  planes  of  the  model  were 
easily  computed  and  drawn  on  the  transparent  coordinate  frame  This 
permitted  easy  visualization  of  the  eigenplanes  and  their  orientation 
with  respect  to  each  other,  to  the  eigenvectors,  and  to  the  eigencone. 


The  first  investigation  of 
the  eigenplane  location  was 
made  by  keeping  the  damping 
ratio,  £,  of  the  system 
constant  while  varying  the 
natural  frequency  of  the 
system,  w^ ,  keeping  the 
real  root  of  the  system 
constant  also.  This  is 
equivalent  to  the  root 
movement  shown  to  the  right 
in  Fig.  7. 


Fig  7  Root  Movement; 

w  variable,  £  constant 
n 


The  root  locations  investigated  were: 
r i  =»  3.0  *  constant  £  *  .707  -  constant 


2,3  ’ 

.5  +  j  . 

.5 

w  *  o  70  7 
n 

(Case 

A) 

2,3  = 

1.0  +  j 

1.0 

w  =*  1 . 414 
n 

(Case 

B) 

=a 

2,3 

2.0  +  j 

2.0 

w  -  2.828 
n 

(Case 

C) 

2,3 

3.0  +  j 

3.0 

w  *  4.242 
n 

(Case 

D) 

Models  of  the  eigenplanes  of  these  cases  were  built  and  from 
these  it  was  seen  that  with  £  constant,  varying  w^  causes  the  eigenplane 
to  '‘twist11  around  the  eigencone.  This  twist  can  be  easily  seen  in  Fig-  8 
which  shows  the  eigenplanes  for  Cases  A,  B,  C  and  D  respectively. 

As  can  also  be  seen  from  Fig.  9,  views  looking  down  the  edges  of  the 

eigenplanes  of  Cases  A,  B,  and  C,  the  angle  that  the  eigenplane  makes  with 

the  eigencone  appears  to  remain  constant  as  w  is  varied.  No  attempt  was 

n 

made  to  analytically  prove  this  angle  to  be  constant,  nor  was  an  attempt 
made  to  derive  a  relationship  of  this  angle  to  the  value  of  £  chosen 
in  the  investigation. 

In  the  next  phase  of  the  investigation,  the  real  root  and  the  real 
part  of  the  complex  roots  were  kept  constant  while  the  imaginary  part  of 
the  complex  roots  was  varied. 
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This  movement  of  the 
roots  is  shown  at  right  in 
Fig.  10. 

The  root  locations  in¬ 
vestigated  were: 

r^  3  3,0  *  constant 


r23  =  1.0+  J  1  0 

r2  3  -  1.0  +  j  2.0 

r2  3  -  1.0  +  j  3.0 

r0  =  1,0  +  j  4.0 

2,3  “  J 


(Case  B ) 
(Case  F) 
(Case  G) 
(Case  H) 


Fig.  10  Root  Movement; 
lm(x  +  jy)  variable, 

Re(x  +  jy)  constant 


In  determining  the  equations  of  the  traces  of  these  complex 
eigenplanes  in  the  coordinate  planes,  it  was  found  that  the  trace  of 
all  complex  eigenplanes,  where  they  intersected  the  E/E  plane,  was 
constant.  This  is  equivalent,  in  three  dimensions,  to  a  rotation  of 
the  complex  eigenplane  about  a  line  in  the  E/E  plane.  This  rotation 
of  the  complex  eigenplane  with  the  variation  of  the  imaginary  part  of 
the  complex  roots  is  shown  pictorially  in  Fig.  11  for  Cases  F,  G  and  H. 

Without  the  use  of  the  three-dimensional  models,  it  can  be  seen  that 
the  variation  of  only  the  imaginary  part  of  the  complex  roots  will  have 
no  effect  on  the  trace  of  the  eigenplane  in  the  E/E  plane.  The  equation 
of  this  trace  is: 

E  +  (r2  +  r3)  E  -  0 

and,  with  r^  and  r ^  complex  conjugates,  the  coefficient  of  the  E  term 
in  the  above  expression  depends  only  on  the  magnitude  of  the  real  part 
of  the  complex  roots. 

The  final  phase  of  this  investigation  of  eigenplane  geometry 
consisted  of  varying  the  system  damping  ratio,  C>  while  holding  the 
natural  frequency  of  the  system,  w  ,  constant. 
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This  movement  of  the  roots  is 
shown  at  the  right  in  Fig.  12. 
The  root  locations  investigated 
were : 

r  =  3.0  *  constant 

w  »  2.83  “  constant 
n 


Fig.  12  Root  Movement; 

£  variable,  w  constant 


r2 , 3 

*  1  414 

+  j2 .45 

£  3  -5 

(Case 

i) 

r2,3 

»  2.0 

+  j2.0 

C  -  .707 

(Case 

C) 

r2 , 3 

-  2.45 

+  jl. 414 

C  *  .866 

(Case 

M) 

r2,3 

-  2.735 

+  j  •  732 

C  *  .966 

(Case 

J) 

Here  again,  in  determining  the  equations  of  the  traces  of  these 

various  complex  eigenplanes,  it  was  found  that  the  trace  of  all  complex 

eigenplanes,  where  they  intersected  one  of  the  coordinate  planes,  was 

constant.  In  this  instance,  the  constant  trace  was  located  in  the 

E/E  plane.  This,  likewise,  is  the  same  as  a  rotation  of  the  complex 

eigenplane  about  a  line  in  the  E/E  plane.  This  rotation  of  the  complex 

eigenplane  as  the  system  C  varies,  with  w  held  constant,  is  shown  in 

n 

Fig.  13,  for  Cases  C,  I,  and  J,  respec ti vel y . 

Further  analysis  of  the  last  two  phases  of  this  investigation  of 
the  geometry  of  the  eigenplane  has  led  to  a  generalization  regarding 
the  orientation  of  the  complex  eigenplane  as  the  location  of  the  complex 
roots  of  a  third  order  system  is  varied  That  is,  as  the  imaginary  part 
of  the  complex  roots  is  varied  from  large  values  to  small  values,  (this 
is  equivalent  to  increasing  the  system  O  the  complex  eigenplane  rotates 
toward  the  surface  of  the  eigencone  and  as  Q  reaches  1.0,  the  complex 
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Case  J  K 

Fig.  13  i 

E.g*nplane  Movement 
7  Variable. ,  u>n  Constant 


2  G 


eigenplane  becomes  tangent  to  the  eigencone  at  the  eigenvector  corres¬ 
ponding  to  the  magnitude  of  these  two  roots  as  they  enter  the  axis  of 
reals.  At  the  moment  of  entry,  the  roots  are  real  and  repeated,  and 
are  associated  with  a  common  eigenvector.  This  may  also  be  shown 
analytically  by  obtaining  the  equation  of  the  eigenplane  tangent  to  the 
eigencone . 

The  general  form  of  a  plane  tangent  to  a  conic  surface  at  a  point  is 


*  ••  •  ^  •• 

where  jF^E,  E,  E)  =  E  -  EE  (the  equation  of  the  eigencone). 

•  * 

Performing  the  partial  dif  ferentiation  at  the  point  E  =  E^,  E  *  E^ , 

E  =»  E  yields 
o 

-E  (E  -  E  )  +  2E  (E  -  E  )  -E  (E  -  E  )  =  0 
o  o  o  o  o  o 

.  «•  /  K  \ 

Let  (Eq,  Eq ,  Eq)  58  —  ,  K,  -r^Kj  ,  any  point  on  an  eigenvector. 

Thus 

r-K(e+7 )  +  ZK(e-k)  f  ^-(E  +  rLK)=  0 
r-  KE  +  Kz  +  2  K  £  - 2  Kz  +  77-  E  +  K '  =  0 

L  L 

Combining  terms  and  dividing  through  by  K  gives: 

•  * 

r.  E+2E+—  =>0 
i  r . 

i 

or 

2 

E  +  2r.  E  +  r.  E  *  0 
l  i 

which  is  the  equation  of  an  eigenplane  defined  by  any  two  repeated 
real  roots ,  r . . 

l 
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A  special  case  of  tangency  between  eigenplane  and  eigencone  exists 
for  the  case  of  three  repeated  real  roots.  Here,  the  line  of  tangency 
corresponds  to  a  single  eigenvector  which  is  the  only  eigenvector  of  the 
system,  and  is  the  eigenvector  for  all  three  roots. 

Extending  this  reasoning  further  helps  explain  the  existance  of  the 
previously  mentioned  pyramid  of  real  root  eigenplanes.  If  a  system 
gain  is  high  enough  so  that  the  root  locus  consists  of  one  real  root 
and  a  complex  pair,  there  will  be  one  real  root  eigenvector  and  a  complex 
eigenplane  in  error  space.  As  the  gain  is  lowered  such  that  the  complex 
roots  just  enter  the  axis  of  reals,  the  complex  eigenplane  is  tangent 
to  the  eigencone  at  the  eigenvector  corresponding  to  the  two  repeated 
real  roots.  As  the  system  gain  is  lowered  further,  the  entering  roots 
will  diverge  on  the  axis  of  reals  and  the  complex  eigenplane  breaks 
into  two  real  eigenplanes,  now  passing  through  the  eigencone  and  inter¬ 
secting  it  at  the  real  root  eigenvectors.  These  two  real  root  eigen¬ 
vectors  now  form  the  pyramid  with  the  original  real  root  eigenvector. 


Method  of  Determining  Roots  to  Locate  an  Eigenplane  in  a  Desired  Location 

A  method  of  fixing  the  eigenplane  in  error  space  has  been  developed 
whereby,  given  any  two  traces  of  the  eigenplane  on  the  coordinate  planes, 
the  required  system  roots  can  be  readily  determined.  As  has  been  pre¬ 
viously  shown,  only  the  complex  conjugate  roots  of  the  system  are  necessary 
to  define  the  complex  eigenplane,  and  the  real  root  of  the  system  can  be 
of  any  magnitude.  The  real  roots,  taken  two  at  a  time,  will  define  the 
eigenplanes  of  a  system  with  three  real  roots. 


The  method  will  first  be  developed  in  general  terms,  and  then  a 
numerical  example  will  be  given.  First,  consider  the  differential 
equation  of  a  third  order  system: 


E*  +  aE  +  bE  +  cE  =»  0 

where:  a  »  r^  +  r^  +  r^ 


b 


rlr2  +  r2r3  +  rlr3 


(111*3) 


c 


r  r  r 
12  3 
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Let  r  *  z 


r2  =  x  +  jy 

r3  *  *2*  *  X  " 

Thus  : 

2  2 

r2  +  r3  3  2x  and  r2r3  ~  x  +  y 

The  equation  of  a  typical  eigenplane  is  given  by: 

M'l  53  E  +  (r2  +  r3)  E  +  r2r3  E  -  0 

Substitution  of  the  above  expressions  for  r^  +  r^  and  r ^r ^  into  yields: 
38  E  +  (2x)  E  +  (x2  +  y2)  E  =  0 

This  expression  will  give  the  following  traces  of  }Ji  in  the  three 
coordinate  planes : 

(In  E/  E  plane)  2xE  +  (x2  4*  y2)  E  3  0  (III-4a) 

(In  E/E  plane)  E  +  2xE  =  0  (III-4b) 

•  •  *  2  2 

(In  E/E  plane)  E  +  (x  +y)E:s0  (III-4c) 

Rewriting  the  above  trace  equations  in  terms  of  the  angle  at  the  origin 
yields : 


X2  4  y2 


(III -5a) 


=  “tan  (3  - 


(III-5b) 


x4  -t-  y2 


(III-5c) 
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Where  a,  3,  and  y  are  defined  in  Fig.  14 


Fig*  14.  Angles  a,  3,  and  y  on  the  coordinate  planes. 

By  determining  any  two  of  the  three  angles,  O',  3,  and  y,  and 
substituting  these  angles  into  the  corresponding  tangent  equations (III-5) , 
a  system  of  two  equations  in  the  two  unknowns,  x  and  y,  are  established. 

On  solving  for  x  and  y,  r^  and  r^  are  determined.  If  y  is  determined  to 
be  complex,  substitution  into  the  expression  r^  ^  s  x  +  jy  yields  r^ 
and  r ^  as  real  roots.  If  y  is  real,  then  r^  and  r^  are  complex  conjugate 
roots . 


Numerical  Example  for  Complex  Eigenplane 

Assume  that  a  complex  eigenplane  is  required  to  be 
that  its  trace  in  the  e/E  plane  is  defined  by  ot  =»  45° 
E/E  plane  by  3  3  60° 

Substituting  these  into  Equations  (IIX-5a)  and  (III-5b) 


fan  4-5 


*  +  3 


^  1 


located  such 
and  in  the 

yields : 


(III*  6a) 


fan  GO 


-L.*  f3~ 


(III-6b) 


Solution  of  (IIl-6b)  directly  determines  x; 

*  _  ,J__  =  .288 

*  "  2VT 
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and  this,  when  substituted  into  (IlT-6a)  yields: 

(+.288)2  +  y2  =  +2  (+.288) 

y?  -  576  -.0833  =  .4927 

°r  y  ■  ^ .4927  =  .701 
Therefore  the  required  roots  r^  and  r^  are: 

r2  3  *  *288  -  J  -701 

Numerical  Example  for  Real  Root  Eigenplane 

Assume  that  an  eigenplane  is  required  to  be  located  such  that  its 
trace  in  the  E/E  plane  is  defined  by  a  *  15  ,  and  in  the  E/E  plane  by 
Y  =  45  . 

Substituting  these  values  into  Equations  (III-5a)  and  (III-5c)  yields: 

tan  IS’--  --  MB 
tan  **+/  =  1 


Substitution  into  the  upper  equation  yields  x  directly: 

_L  -  .  268 

zx 

I  _  .  t  _  I.8GS 

A  '  2  (.26 8)  .S3Q 

(/.8  G5)Z+  ^  =  I 


^  =  I  -3A7S 


3  - 


-ZA75 


±j LSI 3 
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Note 
of  y 


*  +  j  a 

^  I.8G5  +  j  ( j  /.  5  73) 


=  .292 

r3  =  1.865  -j(j,-5‘73J 

r  3.43S 


that  since  and  are  complex  conjugate  roots,  only  one  value 
,  either  positive  or  negative,  need  be  substituted  into  r9  «  =*  x  +  jy. 
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Chapter  IV 


PHASE  TRAJECTORIES  OF  THIRD  ORDER  LINEAR  SYSTEMS 


General 

It  is  well  known  that  the  time  solution  to  the  differential 
equation  of  a  linear  third  order  feedback  control  system  can  be  repre¬ 
sented  by  a  continuous  curve  in  three-dimensional  error  space.  If  the 
differential  equation  of  a  third  order  system, 

E*  +  aE  +  bE  +  cE  -  0  (IV- 1) 

is  solved  by  Laplace  transform  methods,  the  time  solution  obtained  is: 


r/.v  _  £0+^  +  r3)Eo+^ t_0  -r,+  L0  +(V>  QL0  +  (r,  ra) 

t(i)  -  — : - — - £  +  - : - — - 


Eo 


n)('A-r, 


(<V  ri)(S-ri) 


£p  +  (r,f  ra )  £0  +  (r,  rj  £u  -r-,t 
( ri  -  ri)  ( r2  -  rb ) 


(IV-2) 


a  a  • 

where  E  ,  E  ,  and  E  are  the  initial  conditions  of  error,  error  velocity, 
o  o  o 

and  error  acceleration  respectively.  Successive  differentiation  of 
Equation  (IV-2)  with  respect  to  time  yields  similar  expressions  for  E(t) 
and  E ( t) . 

#  •  • 

Substituting  the  values  of  initial  conditions  E  ,  E  ,  and  E  ,  system 

o  o  o 

root  values,  and  time  into  Equation  (IV-2)  and  its  derivatives  yields  the 
coordinates  of  points  which  describe  the  ,,sta.ten,  i.e.,  the  error,  error 
velocity,  and  error  acceleration  of  the  system  at  that  particular  time. 
The  locus  of  all  these  points  is  the  phase  trajectory  of  the  system. 

The  calculation  of  trajectories  by  direct  substitution  into 
Equation  (IV-2)  and  its  derivatives  would  be  a  laborious  process.  This 
simple  but  time  consuming  method  of  solving  for  phase  trajectories 
suggests  solution  by  means  of  electronic  computers. 
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In  the  initial  stages  of  this  investigation,  solution  of  the 
differential  equation  (IV-1)  was  accomplished  utilizing  the  Control  Data 
Corporation  1604  high  speed  solid  state  general  purpose  digital  computer. 
The  digital  computer  program  was  written  in  FORTRAN  language  and  is 
given  in  Appendix  A  along  with  its  flow  diagram.  This  FORTRAN  program 
is  based  on  the  Runge-Kutta  method  of  numerical  integration.  In  order 
to  obtain  the  desired  accuracy  of  solution,  state  points  were  computed 
at  intervals  of  0.01  sec.  but  were  printed  out  at  intervals  of  0.05  sec. 

A  sample  print  out  of  the  digital  solution  is  also  given  in  Appendix  A. 
The  program  was  allowed  to  run  for  a  total  problem  time  of  10-0  sec., 
and  in  practically  all  cases,  at  the  end  of  this  time,  the  values  of 
E,  E,  and  E  were  small  enough  to  consider  that  the  trajectory  had  reached 
the  origin  of  error  space.  In  the  beginning  of  the  investigation,  it  had 
been  anticipated  that  a  digital  to  x-y  plotter  would  be  used  to  obtain 
the  trajectories  in  graphical  form.  However,  this  plotting  equipment  did 
not  become  available,  and  it  was  necessary  to  hand  plot  all  phase 
trajectories . 

In  the  latter  stages  of  the  investigation,  in  order  to  obtain 
trajectories  more  rapidly,  a  Donner  3100  analog  computer  was  used  in 
conjunction  with  an  x-y  servo  plotter.  The  flow  chart  and  scaling  data 
for  the  analog  computer  solution  are  given  in  Appendix  B. 

In  both  the  digital  and  analog  solutions,  plots  of  the  trajectories 
in  the  three  coordinate  planes,  E/E,  E/E,  and  E/E  were  made  and  then 
transferred  to  the  corresponding  planes  of  the  plexiglass  models.  Using 
these  projections,  the  three-dimensional  models  of  the  phase  trajectories 
were  constructed.  Samples  of  the  three  projections  of  a  typical  phase 
trajectory  of  a  three  real  root  system,  as  obtained  from  the  analog 
computer,  are  given  in  Figs.  15  through  17. 

Effect  of  Root  Location  on  Phase  Trajectories 

For  a  system  defined  by  a  given  set  of  roots,  the  phase  trajectories 
in  response  to  various  initial  conditions  are  constrained  by  the  eigen¬ 
vectors  and  eigenplanes  associated  with  the  defining  roots.  Since  it  has 
been  shown  previously  that  eigenvector  and  eigenplane  orientation  are 
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dependent  only  on  root  location,  any  change  in  root  location  will  cause 
a  reorientation  of  eigenvectors  and  eigenplanes  and  so  define  a  different 
system.  The  phase  trajectories  of  this  system  will  be  altered  from 
that  of  the  original  system  in  a  manner  dictated  by  the  new  orientation. 

For  a  system  defined  by  three  real  roots,  it  has  been  previously 
shown  that  any  change  in  root  location  will  effect  the  orientation  and 
size  of  the  pyramid  bounded  by  the  three  eigenplanes-  The  phase  trajec¬ 
tories  will  exhibit  the  same  general  characteristics  in  response  to  initial 
conditions  as  discussed  in  the  following  sections  but  will  be  dependent 
on  the  particular  orientation  and  size  of  the  pyramid 

As  was  previously  indicated,  a  system  defined  by  one  real  root  and 
two  complex  conjugate  roots  has  only  one  eigenvector  and  one  eigenplane. 
Changing  the  location  of  the  real  root  changes  the  time  response  of  the 
system  by  moving  the  eigenvector  along  the  surface  of  the  eigencone 
towards  either  a  slower  or  faster  trajectory.  When  the  complex  roots 
are  changed,  the  eigenplane  is  reoriented  and  the  phase  trajectories  are 
changed  with  regard  to  their  E,  E,  and  E  overshoots.  In  addition  to 
eigenplane  reorientation,  if  the  complex  roots  are  changed  so  as  to 
produce  a  change  in  damping  ratio,  £,  the  general  oscillatory  nature  of 
the  trajectories  is  increased  or  decreased. 

Effect  of  Initial  Conditions  on  Phase  Trajectories 

One  of  the  primary  purposes  of  this  investigation  was  to  study  the 
effect  of  varying  initial  conditions  on  the  phase  trajectories  of  third 
order  linear  systems.  In  particular,  three  general  areas  for  the 
location  of  these  initial  conditions  were  selected- 

1)  Initial  conditions  on  an  eigenvector. 

2)  Initial  conditions  in  an  eigenplane. 

3)  Initial  conditions  not  in  an  eigenplane. 

Each  of  these  will  be  discussed  in  the  following  paragraphs. 

1.  Initial  Conditions  on  an  Eigenvector. 

If  initial  conditions  corresponding  to  any  point  located  on  an  eigen¬ 
vector  are  substituted  into  Equation  (IV-2) ,  two  of  the  three  terms  on 
the  right  side  of  the  equation  vanish.  As  an  example,  let: 
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E  =  —  ,  E  =*  -K ,  and  E  =  Kr.  . 
o  rx  o  o  1 


Substitution  into  Equation  (IV-2)  yields: 


E  (t,  ^  KS^2i3E  £-r* 

Cl'  x  K-r.)(r.- -r.)  ^ 


This  further  simplifies  to 


Lit)  = 


-ci 


From  which 


'  <ancl  Eft)  Kr,  £ 

Hence  it  can  be  seen  that  at  any  time,  t,  the  state  point  of  the 
system  has  the  same  coordinates,  |  “  ,  -K,  Kr^ j  ,  as  any  arbitrary  point 

on  an  eigenvector  and  the  trajectory  is  a  straight  line  coinciding  with 
the  eigenvector.  Figures  18  through  20  are  typical  analog  solutions  of 
a  trajectory  with  initial  conditions  on  a  system  eigenvector. 

This  statement  that  a  phase  trajectory  having  initial  conditions  on 
an  eigenvector  will  remain  on  the  eigenvector  holds  for  both  a  system 
with  three  real  roots,  and  a  system  with  one  real  root  and  two  complex 
conjugate  roots.  This  is  apparent  from  Equation  (IV-2)  which  holds  for 
both  systems. 

As  an  interesting  sidelight  to  this  property  of  phase  trajectories 
with  initial  conditions  on  eigenvectors,  an  unstable  system  with  complex 
roots  in  the  right  half  plane  and  its  real  root  in  the  left  half  plane 
was  simulated  on  the  analog  computer.  This  system  was  given  a  set  of 
initial  conditions  on  an  eigenvector,  and  as  the  solution  was  generated, 

-39- 


-40- 


\[ 

I 


+  E 


4V 


— i - 


_c- 


-- f 


+- 


+  t 


Initial  Conditions 
V  (0.35, -\.0, -3.0)  -  - 

On  Fast  Eigenvector 


\i 


Eg. 18 


E/E  Projection  of  Phase  Trajectory 
.of  System  wi.tllThree  Real  Roots 
T.q=. 0.5  '^.0,  r;  =  3.0 


-E 


i 


i  - 


_  »vMno.tnv  no  jin  nn.i  ouooi  h  on  - 

JHV)  VN  JC1V.V.I  <  n  ' 

*N¥<IWO:i  A11VOW  1  J  \ 


I 


-42- 


the  phase  trajectory  coincided  with  the  eigenvector  until  just  before 
it  reached  the  origin.  Here,  the  instability  of  the  system  caused  the 
development  of  diverging  oscillations  about  the  origin.  Figure  21  is 
the  projection  of  this  trajectory  on  the  E/E  plane. 

2.  Initial  Conditions  in  an  Eigenplane* 

If  initial  conditions  corresponding  to  any  point  in  an  eigenplane 
or  a  complex  eigenplane  are  substituted  into  Equation  (IV-2),  Han  has 
shown  (Ref.  1,  pp.  26-36)  that  the  resulting  phase  trajectory  will 
always  remain  in  the  plane. 

A  phase  trajectory  starting  in  a  complex  eigenplane  will  be  an 
oscillatory  curve  that  spirals  into  the  origin  of  error  space,  always 
remaining  in  the  complex  eigenplane.  Figure  22  shows  the  general 
oscillatory  nature  of  the  trajectories  of  four  different  systems  having 
one  real  root  and  two  complex  conjugate  roots.  The  systens  illustrated 
are : 


Case 

L 

ri 

*  0.5 

r2,3 

*  1  +  j3 

(C  * 

.3) 

Case 

B 

ri 

=“  3.0 

r2,3 

-  1  +  jl 

(C  “ 

.707) 

Case 

C 

ri 

*  3.0 

r2,3 

-  2  +  j2 

(C  * 

.707) 

Case 

A 

ri 

=*  3.0 

r2 , 3 

-  .  5  +  j  .  5 

a  = 

.707) 

Figure  23  is  a  view  of  the  same  model,  but  looking  down  the  edge  of  the 
complex  eigenplane  of  Case  A  Figure  24  shows  the  phase  trajectories  of 
Case  L  only,  looking  into  octant  1  Figure  25  is  a  view  of  the  same 
model  of  Case  L  but  looking  down  the  edge  of  the  complex  eigenplane.  The 
dashed  line  in  Figs.  24  and  25  represents  the  eigenvector  corresponding 
to  the  real  root,  r_(  =*  0.5,  of  Case  L 

In  a  system  having  three  real  and  distinct  roots  the  phase  trajec¬ 
tories  will  also  remain  in  an  eigenplane  if  the  initial  conditions  are 
located  in  this  eigenplane.  However,  the  trajectories  will  not  have  the 
oscillatory  nature  as  in  Figs.  22  through  25,  but  will  be  curves  that 
become  asymptotic  to  the  slower  of  the  two  eigenvectors  which  define 
the  eigenplane. 
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Figure  26  shows  a  model  constructed  of  a  number  of  phase  trajec¬ 
tories  of  a  system  with  three  real  roots  ,  i  ^  -  0  5 ,  ~  1  0 , 

r^  =*  2.0,  all  having  initial  conditions  in  the  Third  eigenplane,  defined 
by  the  slow  and  intermediate  eigenvectors.  Note  here  that  two  of  the 
phase  trajectories,  (A)  and(B),  are  straight  lines  in  octants  4  and  21, 
indicating  that  the  initial  conditions  for  these  trajectories  were 
chosen  on  the  slow  eigenvector  This  figure  readily  illustrates  the 
above  statement  that  the  phase  trajectories  become  asymptotic  to  the 
slower  of  the  two  defining  eigenvectors. 

Figure  27  is  another  view  of  the  same  model  but  looking  down  the 
edge  of  the  Third  eigenplane,  and  showing  that  the  phase  trajectories 
all  remain  in  the  eigenplane. 

Figure  28  shows  phase  trajectories  of  a  system  having  three  real 
roots,  r^  *  0.333,  r^  35  1.0,  r.^  35  3.0.  In  this  figure,  Trajectory  (A) 
has  initial  conditions  in  the  First  eigenplane  and  becomes  asymptotic 
to  the  intermediate  eigenvector,  Trajectory  (B)  has  initial  conditions 
in  the  Second  eigenplane  and  becomes  asymptotic  to  the  slow  eigenvector, 
and  Trajectory  (C)  has  initial  conditions  in  the  Third  eigenplane  and 
also  becomes  asymptotic  to  the  slow  eigenvector. 

3.  Initial  Conditions  not  in  an  Eigenplane. 

Han  and  Thaler  have  stated  (Ref,  3)  that  eigenplanes  "subdivide 
the  (error)  space  into  regions  so  that  no  phase  trajectories  can  pass  out 
of  one  region  and  into  another.  Thus  the  (eigenplanes)  act  as  boundaries 
which  funnel  the  phase  trajectories  into  the  origin". 

For  a  system  having  three  real  and  distinct  roots,  this  statement 
can  be  modified  to  read  that  "the  pyramid  of  eigenplanes  subdivides 
error  space  into  regions  so  that  phase  trajectories  starting  inside 
the  pyramid  will  remain  inside  and  those  phase  trajectories  originating 
outside  the  pyramid  will  never  enter  the  pyramid  of  eigenplanes.  Thus 
the  pyramid  acts  as  a  boundary  which  funnels  the  phase  trajectories  into 
the  origin". 

It  is  important  to  note  that  not  only  does  the  pyramid  funnel  the 
phase  trajectories  into  the  origin,  but  it  funnels  them  in  asymptotic 
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to  the  slow  eigenvector,  no  matter  how  close  the  initial  conditions 
might  be  to  the  fast  or  intermediate  eigenvectors  Figure  29  shows 
two  phase  trajectories  of  a  three  real  root  system  with  roots: 
r^  =*  0.333,  r^  ~  1.0,  r^  »  3.0,  whose  initial  conditions  are  inside 
the  pyramid  of  eigenplanes.  It  can  readily  be  seen  that  both  trajec¬ 
tories  approach  the  slow  eigenvector  asymptotically  even  though  (A)  has 
its  initial  conditions  near  the  intermediate  eigenvector  and  (B)  has 
its  initial  conditions  approximately  halfway  between  the  intermediate  and 
fast  eigenvectors. 

Figure  30  is  an  illustration  of  a  number  of  phase  trajectories 
whose  initial  conditions  are  all  outside  the  pyramid  of  eigenplanes. 

Here  again  it  can  be  seen  that  all  phase  trajectories  originating 
outside  the  eigenplanes  approach  the  origin  of  error  space  asymptotic  to 
the  slow  eigenvector.  Three  of  these  trajectories  exhibit  interesting 
properties.  Trajectory  (A)  originates  in  octant  4  near  the  fast  eigen¬ 
vector,  but  does  not  approach  the  origin  along  the  slow  eigenvector  in 
the  same  octant.  Instead  it  penetrates  into  octants  3  and  2,  then 
finally  approaches  the  origin  asymptotic  to  the  slow  eigenvector  in 
octant  2'.  Trajectory  (B)  originates  in  octant  4',  but  instead  of  going 
to  the  slow  eigenvector  in  the  adjacent  octant  4,  it  journeys  through 
octants  4,  3,  2,  and  finally  approaches  the  origin  along  the  slow  eigen¬ 
vector  in  octant  2'.  Trajectory  (C)  originates  near  the  intermediate 
eigenvector  in  octant  4,  approaches  this  intermediate  eigenvector, 
parallels  it  for  a  short  time,  and  then  finally  becomes  asymptotic  to 
the  slow  eigenvector  very  near  the  origin  in  octant  4. 

These  unusual  excursions  through  error  space  of  trajectories  whose 
initial  conditions  are  outside  the  pyramid  of  eigenplanes  are  results 
of  the  basic  property  of  eigenplanes,  i.e.,  eigenplanes  act  as  boundaries 
on  phase  trajectories.  When  the  initial  conditions  are  between  any  two 
eigenplanes,  the  phase  trajectories  must  always  remain  between  them  until 
the  trajectories  reach  the  slow  eigenvector  which  exists  between  or  on 
these  bounding  eigenplanes. 
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In  the  investigation  of  the  effect  of  initial  conditions  not  in  an 
eigenplane,  more  attention  was  focused  on  those  systems  having  one  real 
root  and  two  complex  conjugate  roots.  A  number  of  interesting  results 
were  obtained 

The  first  complex  root  system  which  was  studied  had  roots  of: 
r^  3  3.0,  ^  9  *25  jr  j5.0.  This  system  was  very  lightly  damped, 

having  a  value  for  £  of  0.05.  As  might  be  expected,  the  phase  trajec¬ 
tory  was  highly  oscillatory.  This  can  readily  be  seen  in  Fig  31. 

In  this  particular  three-dimensional  model,  the  scaling  of  E,  E,  and  E 
values  of  points  on  the  phase  trajectory  were  modified  in  order  to 
obtain  a  model  capable  of  analysis.  Here,  for  one  coordinate  axis  unit, 

E  38  2 ,  E  =  4,  and  E  *  1.  When  a  random  initial  condition,  Point  (A), 
not  in  the  complex  eigenplane,  was  applied  to  this  system,  the  phase 
trajectory  very  quickly  became  asymptotic  to  the  complex  eigenplane. 

Figure  32,  a  view  looking  down  the  E  axis,  shows  that  the  phase  trajec¬ 
tory  has  become  asymptotic  to  the  complex  eigenplane  after  less  than  one 
oscillation.  Here  it  should  be  noted  that  this  model  pictured  in  Figs.  31 
and  32  does  not  show  the  last  five  oscillations  of  the  complete  trajec¬ 
tory,  and  they  would  continue  from  Point  (B)  and  eventually  reach  the 
origin  of  error  space. 

The  next  complex  root  system  studied  had  roots  at  r^  =05, 
r2  3  *  1.0  +  j3  0  The  system  was  lightly  damped,  having  C  3  0.3. 

The  general  nature  of  this  system's  oscillatory  trajectories  was  seen 
previously  in  Fig.  24  In  this  investigation,  a  variety  of  initial 
conditions  not  in  the  complex  eigenplane  were  applied  and  analyzed 

The  first  initial  condition  applied  was  a  step  input  In  this 
instance,  the  step  had  a  magnitude  of  4  0  radians.  Curve  (A)  of 
Fig.  33  is  the  resulting  phase  trajectory  Also  shown  in  Fig  33  is 
the  real  root  eigenvector  represented  by  the  dashed  line,  and  two 
trajectories  having  initial  conditions  in  the  complex  eigenplane, 

Curves  (B)  and  (B').  Figure  33  indicates  that  for  this  system,  the 
phase  trajectory  starting  outside  the  complex  eigenplane  does  not 
approach  the  origin  of  error  space  asymptotic  to  the  complex  eigenplane, 
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but  instead  appears  to  become  asymptotic  to  the  system  eigenvector 
Figure  34  is  another  view  of  the  same  model ,  but  looking  up  the  eigen¬ 
vector  in  octant  4~  Note  here  the  similarity  of  the  general  shape  of  the 
phase  trajectories  having  initial  conditions  in  and  out  of  the  complex 
eigenplane. 

The  next  initial  conditions  applied  to  this  system  was  an  initial 
velocity  of  magnitude  2.0  rad/sec  .  Curve  (A)  of  Fig  35  is  the  resulting 
phase  trajectory.  Also  shown  are  the  system  eigenvector  and  the  two 
trajectories  in  the  complex  eigenplane  as  before.  Fig.  36  is  a  view 
of  the  same  model  looking  down  the  edge  of  the  complex  eigenplane. 

Although  it  is  difficult  to  tell  from  Fig,  36,  the  phase  trajectory  did 
not  immediately  become  asymptotic  to  the  complex  eigenplane,  as  it  did 
in  the  preceding  system  having  £  86  0.5.  Figure  37  is  another  view  of 
the  same  model  looking  down  the  eigenvector  in  octant  2*.  Note  again  the 
similar  curvature  of  the  phase  trajectories  having  initial  conditions 
in  and  out  of  the  complex  eigenplane. 

Next,  an  initial  condition  very  near  the  eigenvector  was  applied  to 
the  system.  As  can  be  seen  from  Fig.  38,  the  phase  trajectory  quickly 
became  asymptotic  to  the  eigenvector  instead  of  to  the  complex  eigenplane. 

Finally  an  initial  condition  that  was  equidistant  from  the  eigen¬ 
vector  and  the  complex  eigenplane  was  applied  to  the  system  Figure  39 
shows  the  resulting  phase  trajectory  as  Curve  (A)  and  the  two  trajectories 
in  the  eigenplane,  curves  (B)  and  (B1)  as  well  as  the  system  eigenvector. 
As  can  be  seen  in  Fig,  39,  the  phase  trajectory  again  becomes  asymptotic, 
not  to  the  complex  eigenplane,  but  to  the  eigenvector,  developing  the 
characteristic  spiral  around  the  eigenvector  as  was  seen  previously  in 
Figs.  33  and  38.  Figure  40  is  a  view  of  this  model  looking  up  the 
eigenvector  in  octant  4.  The  similarity  of  the  phase  trajectories 
with  initial  conditions  in  and  out  of  the  complex  eigenplane  noted  here 
and  in  previously  discussed  views  looking  up  the  eigenvector  indicated 
that  there  must  be  some  relation  between  the  phase  trajectories  and  the 
eigenvector.  This  led  to  the  next  phase  of  the  investigation. 
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However,  before  proceeding  further,  the  results  of  this  phase  of 
the  investigation  indicate  that  for  systems  having  one  real  root  and  two 
complex  roots,  the  phase  trajectories  of  the  systems  that  are  very  lightly 
damped  tend  to  become  asymptotic  to  the  complex  eigenplane.  For  systems 
with  moderate  to  heavy  damping,  the  trajectories  become  asymptotic  to 
the  system  eigenvector  as  they  approach  the  origin  of  error  space. 

Interrelation  of  Phase  Trajectories,  Eigenplanes,  and  Eigenvectors 

In  the  study  of  phase  trajectories  of  systems  having  one  real  and 
two  complex  conjugate  roots,  it  was  noted  that  a  striking  similarity 
existed  between  trajectories  with  initial  conditions  inside  and  outside 
the  complex  eigenplane  when  the  trajectories  were  viewed  looking  along 
the  system  eigenvector,  as  in  Fig*  40  Further  investigation  led  to 
the  following  development  of  an  interrelation  among  phase  trajectories, 
eigenplanes,  and  eigenvectors. 

It  can  be  seen  from  the  mathematics  of  Chapters  II  and  III,  and 
also  in  Figs.  24  and  25,  that  the  system  eigenvector  pierces  the  complex 
eigenplane  at  the  origin  of  error  space.  Mathematically,  the  complex 
eigenplane  can  be  translated  along  the  eigenvector,  keeping  the  origin 
of  the  complex  eigenplane  coincident  with  the  eigenvector,  and  keeping 
the  complex  eigenplane  oriented  parallel  to  its  original  position.  A 
set  of  initial  conditions  can  be  applied  to  the  system  corresponding  to 
the  point  on  the  eigenvector  containing  the  origin  of  the  translated 
eigenplane.  If  the  complex  eigenplane  is  then  allowed  to  move  along 
the  eigenvector  at  the  same  velocity  as  the  system  state  point,  this 
state  point  will  always  coincide  with  the  origin  of  the  complex  eigen¬ 
plane  as  these  points  move  toward  the  origin  of  error  space  Extending 
this  further,  let  the  complex  eigenplane  be  translated  again  along  the 
eigenvector,  and  impose  a  set  of  initial  conditions  on  the  system  such 
that  the  initial  state  point  lies  in  the  translated  eigenplane  some 
distance  from  the  translated  originc  If  the  complex  eigenplane  is 
allowed  to  move  toward  the  origin  along  the  eigenvector  as  before,  the 
state  point  now  traces  the  phase  trajectory  in  three-dimensional  error 
space.  What  has  really  happened  is  that  the  state  point  has  never  left 
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the  translated  complex  eigenplane  The  trace  it  has  made  in  the  complex 
eigenplane  is  exactly  the  same  as  the  phase  trajectory  described  by  the 
state  point  had  the  initial  conditions  been  applied  in  the  original 
untranslated  complex  eigenplane  at  the  same  distance  from  the  origin  as 
the  present  initial  conditions  were  from  the  eigenvector  in  the  trans¬ 
lated  plane.  Hence  it  can  be  seen  that  any  phase  trajectory  is  the 
result  of  the  superposition  of  two  separate  sets  of  Initial  conditions, 
one  set  on  the  eigenvector  and  the  other  set  in  the  complex  eigenplane. 

In  a  system  with  three  real  roots,  the  same  results  are  obtained  if 
an  eigenplane  is  translated  along  the  eigenvector  not  defining  the  eigen¬ 
plane  . 

Hence,  any  phase  trajectory  in  three-dimensional  error  space  can 
be  described  as  the  superposition  of  the  two  trajectories  resulting  from 
initial  conditions  on  an  eigenvector,  and  initial  conditions  in  either  a 
complex  or  real  root  eigenplane. 

This  hypothesis  was  then  checked  using  the  analog  computer  set-up 
as  given  in  Appendix  C.  In  Fig  C-l,  the  upper  system  was  given  initial 
conditions  corresponding  to  a  point  lying  on  the  eigenvector,  The  lower, 
or  "primed11  system  was  given  a  set  of  initial  conditions  corresponding  to 
a  point  offset  from  the  eigenvector  by  a  distance  equivalent  to  the  initial 
conditions  in  the  complex  eigenplane.  The  difference  of  the  two  systems 
was  then  plotted  by  x-y  recorder  and  compared  with  the  phase  trajectory 
in  the  complex  eigenplane. 

Figure  41  is  the  E/E  projection  of  the  phase  trajectory  in  the 

•  •  * 

complex  eigenplane  and  Fig  42  is  the  E/E  projection  of  the  phase  tra¬ 
jectory  of  the  "difference"  system  The  validity  of  the  above  hypothesis 
can  be  seen  by  a  comparison  of  the  identical  projections  of  the  phase 
trajectories  of  Figs  -  41  and  42  c 

Complex  Root  Phase  Trajectories  in  Response  to  a  Step  Input 

In  Fig.  33,  the  phase  trajectory  of  a  system  with  one  real  root  and 
a  pair  of  complex  conjugate  roots  is  shown  for  a  step  input  By  noting 
the  tight  spiral  of  the  trajectory  around  the  eigenvector,  it  must  be 
concluded  that  for  this  particular  case,  the  trajectory  intersects  the 
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eigencone.  Since  trajectories  of  this  type  could  prove  of  great  impor¬ 
tance  in  switching  applications,  it  was  decided  to  investigate  additional 
trajectories  with  step  inputs  to  see  if  a  definite  pattern  could  be 
established  as  to  if,  where,  and  when  the  trajectories  intersect  the 
eigencone . 

The  three-dimensional  models  which  had  proved  invaluable  heretofore 
could  not  be  constructed  precisely  enough  tc  provide  anything  more  than 
a  general  description  of  the  trajectory  However,  from  the  models  it  was 
established  that  all  complex  trajectories  do  not  intersect  the  eigencone. 
A  trajectory  of  this  type  is  shown  in  Fig  31.  Also,  a  study  of  Figs.  5 
and  33,  indicates  that  if  a.  trajectory  does  intersect  the  eigencone,  it 
will  occur  at  a  faster  eigenvector  than  that  of  the  system 

After  having  exhausted  the  possibilities  of  the  three-dimensional 

models,  an  analytical  approach  was  adopted  If  the  initial  conditions 

of  a  phase  trajectory  are  chosen  such  that  Eq  «  *  0,  the  initial 

condition  E  is  an  initial  step  and  can  be  considered  as  a  step  input, 
o 

Introducing  these  initial  conditions  into  Equation  (IV-2): 


(IV-3) 
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Taking  the  first  and  second  derivatives  of  the  above  expression* 
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Substituting  the  above  expressions  for  E(tl,  E(t),  and  E(t)  into 
Equation  (11-14),  the  equation  of  the  eigencone* 


=  O 


(IV- 6) 


This  equation  satisfies  the  conditions  that  exist  when  and  if  the  phase 
trajectory  in  response  to  a  step  input  intersects  the  eigencone.  It 
can  be  seen  from  this  equation  that  i f  an  intersection  occurs,  the  time 
of  intersection  is  independent  of  the  size  of  the  step  input.  Further¬ 
more,  on  obtaining  a  time  of  intersection,  the  equation  becomes  linear 
and  indicates  that  all  points  of  intersection  are  along  a  straight  line 
on  the  surface  of  the  eigencone.  This  line,  by  definition,  must  be 
along  the  same  eigenvector 

In  order  to  obtain  a  solution  to  Equation  (IV~6) ,  the  equation  is 
simplified  as  shown  in  Appendix  D  to  obtain  the  transcendental  equation: 


Cos  3  t  +  o  Sin  3  t  56  e 


(1V-7) 


where:  a  =  real  root 


a  3  real  part  of  complex  root 
3  *  imaginary  part  of  complex  root 


•  i 
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If  a  time  of  intersection  exists  other  than  the  trivial  case 
of  t  =*  0,  it  can  be  obtained  from  a  graphical  solution  of  the  above 
equation.  If  a  solution  exists,  in  all  except  special  cases,  it  is 
possible  to  obtain  more  than  one  solution  if  the  plot  is  continued 
beyond  the  first  intersection.  Figure  33  is  a  trajectory  which  would 
produce  multiple  intersections  as  noted  by  the  tight  spiral  around 
the  eigenvector. 
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In  order  to  avoid  plotting  Equation  (IV- 7)  for  systems  with  no 
intersections,  R.  L.  Ashford  (Ref,  6)  has  devised  a  criteria  whereby 
it  is  possible  to  determine  whether  or  not  a  solution  exists.  In 
essence,  it  amounts  to  expanding  both  sides  of  Equation  (IV-7)  in  a 
Maclaurin  series  and  taking  the  coefficient  of  the  first  power  term  in 
order  to  determine  the  slope  of  each  curve.  The  slope  of  the  curve  for 
the  left  side  of  the  equation  is  given  by  a 3  while  that  of  the  right  is 
(a-a).  Ashford  states  that  when  (a- O')  is  less  than  zero,  there  will 
always  be  an  infinite  number  of  solutions,  and  when  (a- a)  is  greater 
than  zero,  a  non-trivial  solution  will  exist  provided  tff3  is  greater 
than  (a-a).  Ashford’s  criteria,  then,  states  that  a  phase  trajectory 
in  response  to  a  step  input  will  intersect  the  eigencone  whenever: 

1.  (a-a)  <  0 

2.  (a-cv)  >  0  and  0“3  >  (a -ck) 

If,  through  Ashford’s  criteria,  an  intersection  is  indicated,  a 
graphical  solution  of  Equation  (IV-7)  will  determine  the  time  of  inter¬ 
section.  By  substituting  this  value  of  time  into  Equations  (IV-3) , 
(IV-4),  and  (IV-5),  the  point  of  intersection  of  the  trajectory  with 
the  eigencone  is  determined  and  hence,  the  eigenvector. 
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Chapter  V 


CONCLUSIONS 

From  the  results  of  this  investigation,  the  following  conclusions 
can  be  made  regarding  the  eigenvector  geometry  and  phase  trajectories  of 
third  order  linear  feedback  control  systems: 

1)  The  orientation  of  system  eigenvectors  in  third  order  error 
space  if  a  function  of  the  location  of  the  roots  of  the  system  on  the 
complex  s-plane. 

2)  The  locus  of  all  eigenvectors  associated  with  a  third  order 
system  is  a  conic  surface  describing  two  right  elliptical  cones  located 
in  octants  1,  2,  3'  and  4*  with  apexes  at  the  origin  of  error  space. 

3)  An  eigenplane  can  be  fixed  at  a  desired  orientation  in  error 
space  by  proper  choice  of  system  roots. 

4)  The  orientation  of  the  complex  eigenplane  in  error  space  is 
completely  independent  of  the  real  root  of  the  system  and  depends  only 
on  the  location  of  the  complex  conjugate  roots  on  the  complex  s-plane. 

5)  The  complex  eigenplane  never  intersects  the  eigencone  except  at 
the  origin,  however,  it  can  become  tangent,  and  its  point  of  tangency  is 
along  the  eigenvector  associated  with  a  pair  of  repeated  real  roots. 

6)  The  three  eigenplanes  of  a  system  having  three  real  and  distinct 
roots  form  a  pyramid  interior  to  the  eigencone.  This  pyramid  serves  to 
funnel  phase  trajectories  into  the  origin  of  error  space. 

7)  Phase  trajectories  of  systems  with  three  real  roots  having 
initial  conditions  in  an  eigenplane  will  remain  in  this  plane  and  will 
approach  the  origin  of  error  space  asymptotic  to  the  slowest  of  the  two 
defining  eigenvectors . 

8)  Phase  trajectories  of  systems  with  one  real  and  two  complex 
conjugate  roots  having  initial  conditions  in  the  complex  eigenplane  will 
remain  in  the  complex  eigenplane. 

9)  Phase  trajectories  of  moderate  and  heavily  damped  complex  root 
systems,  having  initial  conditions  not  in  the  complex  eigenplane,  spiral 
into  the  origin  of  error  space  asymptotic  to  the  system  eigenvector. 
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10)  Phase  trajectories  of  very  lightly  damped  complex  root  systems, 
having  initial  conditions  not  in  the  complex  eigenplane,  will  spiral 
toward  the  origin  asymptotic  to  the  complex  eigenplane. 

11)  Phase  trajectories  of  a  system  with  three  real  roots  whose 
initial  conditions  are  not  in  an  eigenplane  cannot  penetrate  an  eigen¬ 
plane  and  will  approach  the  origin  of  error  space  asymptotic  to  the  slow 
eigenvec  tor . 

12)  Any  third  order  phase  trajectory  in  response  to  any  initial 
conditions  can  be  described  by  the  superposition  of  the  two  phase  trajec¬ 
tories  resulting  from  initial  conditions  on  an  eigenvector  and  initial 
conditions  in  either  a  complex  or  real  root  eigenplane.  (See  Chapter  IV, 
PP-  67  » 68  ) 

13)  If  the  phase  trajectory  of  a  system  with  one  real  and  a  pair 
of  complex  conjugate  roots  intersects  the  eigencone  in  response  to  a 
step  input,  the  time  of  intersection  is  independent  of  the  size  of  step 
and  the  intersection  will  occur  along  a  faster  eigenvector  than  the 
original  system  eigenvector. 
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Chapter  VI 


RECOMMENDATIONS 

From  the  results  of  this  thesis  it  can  be  seen  that  there  are 
certain  areas  of  interest  which  warrant  further  investigation.  A  logical 
extension  of  this  thesis  would  be  an  attempt  to  project  the  results 
obtained  herein  for  third  order  systems  into  the  more  general  case  of 
Nth  order  systems.  The  existence  of  hypersurfaces  in  N-dimensiona 1 
error  space  corresponding  to  the  three-dimensional  eigenplanes  has 
been  well  established  in  the  literature.  While  not  yet  proven,  it 
seems  reasonable  to  assume  that  there  must  also  exist  in  N-dimensions 
a  hypersurface  which  corresponds  to  the  eigencone.  Thus  it  should  be 
possible  to  describe  the  characteristics  of  phase  trajectories  of  Nth 
order  systems  in  terms  of  the  various  hypersurfaces  defined  by  the 
N-dimensional  eigenvectors . 

The  results  of  this  thesis  could  also  have  applications  in  the 
design  of  discontinuous  control  systems.  The  literature  contains  a 
great  deal  of  information  on  various  schemes  which  utilize  hypersurfaces 

as  switching  surfaces,  but  none  of  these  hypersurfaces  are  comparable  to 

* 

the  eigencone.  However,  in  the  N-dimensional  case,  it  would  have  to  be 
proven  mathematically  that  under  certain  conditions  the  phase  trajectories 
intersect  an  N-dimensional  eigencone  and  that  eigencone  switching  could 
be  effected. 

The  effect  of  non-linearities  on  the  phase  trajectories  of  third 

order  and  higher  systems  is  another  field  of  investigation  to  which  the 

results  of  this  thesis  can  also  be  extended.  The  regions  of  three- 

dimensional  error  space  in  which  a  system  is  linear  can  be  defined  by 

various  planes  which  establish  the  boundaries  of  the  linear  region.  For 

example,  in  the  case  of  velocity  saturation,  the  linear  region  is  defined 
* 

by  the  planes  E  ~  +  K,  where  K  is  a  constant.  Similarly,  the  linear 
region  in  the  case  of  torque  saturation  is  established  by  the  planes 
*E  3  +  K.  In  N-dimensions  it  seems  reasonable  to  assume  that  corresponding 
hypersurfaces  exist  which  also  establish  the  boundaries  of  these  linear 
regions,  and  their  effect  on  phase  trajectories  has  not  been  adequately 
described. 
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A?r TT"DIX  A 

DIGITAL  C^:,?UTGR  FLO'"  DIAGRAM,  PPOG^A'  , 
AND  SAMPLE  PRINT  OUT  OP  THIRD  ORDER 
DIFFERENTIAL  EQUATION  SOLUTION 


START 


Fig.  A-l  Flow  Diagram  of  Digital  Computer  Program 
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FORTRAN  PROGRAM  FOR  DIGITAL  COMPUTER  SOLUTION  OF  THIRD  ORDER  DIFFERENTIAL  EQUATION 


C 


6 

1  3 
4 
9 

1  1 

7 


3 


2 

1 

3 


PROGRAM  PHATRA.J 

D I E  T  7  - WE  S  T  8 ROOK  THESIS  COMPLEX  CASE  F 

dimension  rm 

READ  6,  DEL  TAT  »N 
FORMAT  (  F  5 . 2  »  SX,  13) 

READ  4,  T,  ( Y ( I ) »  1=1,3) 

FORMAT  (F7.3,  3 (  8  X  ,  F  8 . 4  )  ) 

11  =  1 

CALL  RKUTTA  (3,  T,  Y,  DELTA!) 

IF  ( I  I -N )  7,  8,  8 

I  1  =  I  1  +  I 
T  =  I  +DELTA  T 

IF  (T-IO.O)  l l ,  ll,  13 
T  = T  +  DELTAT 

WRITE  OUTPUT  TAPE  8,  9,  T,  (Y(I),  1=1,3) 

FORMAT  (  311  T=  F7.3,4X,3HE=  F 8 . 4 , 4X , 4 HE l =  F  8 . 4 , 4  X  ,  4Ht2=  F8.4) 

GO  TO  9 
END 

SUBRGUT I NE  RKUTTA ( NUMBER, X VAR , YVARS , S  TEP ) 

DIMENSION  YVARS ( 30 ) , AK ( 4 , 30 ) , D Y ( 30 ) , YC ( 30 ) , C ( 4 ) 

C  (  I  )  =  0 . 0 

C ( 2  )  =0 . 5 

C  (  3  )=0.‘3 

C (4 )  =  1  .0 

DO  l  1=1,4 

X  C  =  X  V  A  R  *  c (  1  ) -STEP 

DO  2  J= l, NUMBER 

YC ( J )  = Y V  AK  S ( J )  +  C(  I  ) *AK ( 1-1 »  J ) 

CALL  UERIV(XC,YC»DY) 

DO  1  J  =  l , NUMBER 

AK ( I , J )  =  STEP  »  DY ( J ) 

DO  3  J  =  1 , NUMBER 

YVARS! J)  =YVARS(J)+!AK( 1,J)+2«*AK(2»J)+2.*AK(3,J)fAK(4,J) )/6. 

RE  TURN 

END 

SUBROUTINE  DERIV  (T,  Y,  DY ) 

DIMENSION  Y ( 3 ) ,  DY(3) 

A=4  *0 
8=3.8 
C  =  l  .8  - 
DY!  I  )  =  Y ! 2  ) 

DY ( 2 ) =Y ( 3 ) 

DY (  3  )  =-C  * Y (  1  )-H*Y ( 2 )-A*Y ( 3 ) 

RE  TURN 

END 

tND 


Sample  Digital  Print  Out  -  PROGRAM  PHATRAJ 


Initial  Conditions  :  on  Eigenvector 


T= 

.000 

E=* 

-.6667 

E 1  = 

2 . 0000 

E2=> 

'6.0000 

r= 

.050 

fc  = 

-.5733 

El  = 

1.7214 

£2  = 

-5.  1642 

I: 

:t?8 

E* 

-.4939 

4  1 

-. 3659 

b  1  = 
b  t* 

L-iWS 

r^: 

-4.4449 
~3. 8250 

T  = 

.200 

e  = 

£1  = 

1.0976 

£2  = 

-3.2929 

T  = 

.250 

e= 

-.3149 

El  = 

.9447 

£2  = 

-2.834? 

T  = 

.300 

E  = 

-.2711 

E 1  = 

.8131 

F2  = 

-2.4394 

T  - 

.350 

£  = 

-.2333 

£1  = 

.6999 

£2  = 

-2.0996 

T  = 

.400 

£  = 

-.2000 

b  1  = 

.6024 

£2  = 

-1 . 807? 

T  = 

.450 

£  = 

-.  ! 729 

£1  = 

.5  185 

£2  = 

-1 . 5554 

T- 

.500 

6  = 

-.1488 

bl  = 

.4463 

£2  = 

-1.3388 

T  = 

.550 

E  = 

-.1231 

£1  = 

.3841 

£2  = 

-1.1523 

T  = 

.600 

b  = 

-. 1  102 

E  l  = 

.3306 

£2  = 

-.9913 

T  = 

.650 

t  = 

-.0949 

£  l  = 

.2846 

£2  = 

-.8536 

T  = 

.700 

t  = 

-.0817 

ei= 

.2449 

£2  = 

-.7347 

T  = 

.750 

E  = 

-.0703 

b  1  = 

.2  108 

£2  = 

-.6324 

T  = 

.800 

E  = 

-.0605 

61  = 

.  1814 

E2  = 

-.5443 

T  =  • 

.850 

b  = 

-.0521 

£1  = 

.  1562 

L2  = 

-.4685 

T  = 

.900 

b  = 

-.0448 

£1  = 

.  1  344 

1 2  = 

-.403? 

T  = 

.950 

b  = 

-.0386 

£1  = 

.  1  157 

£2  = 

-.3471 

T  =■ 

1.000 

E  = 

-.0332 

£1  = 

.0996 

E2  = 

-.2987 

1  = 

1.050 

b  = 

-.0286 

bl  = 

.0857 

£2  = 

-.2571 

r= 

1.  100 

E  = 

-.0246 

bl  = 

.0736 

£2  = 

-.2213 

T  = 

1.150 

b  = 

-.0212 

E  1  = 

.0635 

£2  = 

-.  1905 

T = 

1.200 

b  = 

-.0182 

bl  = 

.0547 

£2  = 

-.  1639 

T- 

1.250 

E- 

-.0  157 

£1  = 

.0470 

E2  = 

-.1411 

T  = 

1.300 

F= 

—  . 0  13  5 

E  1  = 

.0405 

£2  = 

-.1214 

T  = 

1.350 

b  = 

- .0  1  16 

£1  = 

.0349 

£2  = 

-.1045 

T  = 

1.400 

b  = 

-.0100 

£1  = 

.0300 

£2  = 

-.0900 

T  = 

1.450 

6= 

-.0086 

b  1  = 

.0258 

E2  = 

-.0774 

T= 

1.500 

E- 

-.0074 

£1  = 

.0222 

£2  = 

-.0667 

T  = 

1.550 

£  = 

-.0064 

£1  = 

.0191 

E2  = 

-.0574 

T  = 

1.600 

E  = 

-.0055 

E  l  = 

.0165 

£2  = 

-.0494 

T  = 

1 .650 

E  = 

-.0047 

£1  = 

.0142 

£2  = 

-.0425 

T = 

1.700 

£  = 

-.004  l 

E  1  = 

.0122 

£2  = 

-.0366 

T  = 

1.750 

£'  = 

-.0035 

£1  = 

.0105 

E2  = 

-.0315 

T = 

1.800 

£  = 

-.0030 

E  1  = 

.0090 

62  = 

-.0271 

T- 

1.850 

E  = 

-.0026 

£1  = 

.0078 

E2  = 

-.0233 

T  = 

1.900 

b  = 

-.0023 

El  = 

.0067 

£2  = 

-.0201 

T = 

1.950 

b  = 

-.0019 

£1  = 

.0058 

E2  = 

-.0173 

T = 

2.000 

b  = 

-.0017 

E 1  = 

.0050 

£2  = 

-.0149 

T  = 

2.050 

£  = 

-.0014 

E  1  = 

.0043 

c2  = 

-.0128 

T  = 

2.100 

E= 

-.0012 

E 1  = 

.0037 

£2  = 

-.0110 

T  = 

2.150 

E  = 

-.0011 

£1  = 

.0032 

E2  = 

-.0095 

T  = 

2.200 

£  = 

-.0009 

£1  = 

.0027 

£2  = 

-.0082 

T  = 

2.250 

E  = 

-.0008 

E 1  = 

.0024 

£2  = 

-.0070 

T  = 

2.300 

£  = 

-.0007 

E  1  = 

.0020 

£2  = 

-.0060 

T  = 

2.3  50 

b  = 

-.0006 

£1  = 

.0017 

E2  = 

-.0052 

T  = 

2.400 

E  = 

-.0005 

E 1  = 

.0015 

£2  = 

-.0045 

T= 

2.450 

b  = 

-.0004 

£1  = 

.0013 

£2  = 

-.0039 

T  - 

2.500 

£  = 

-.0004 

£1  = 

.0011 

£2  = 

-.0033 

T  = 

2.550 

fc  = 

-.0003 

£1  = 

.0010 

£2  = 

-.0029 

r  = 

2.600 

E- 

-.0003 

£1  = 

.0008 

£2  = 

-.0025 

T= 

2.650 

E- 

-.0002 

£1  = 

.0007 

£2  = 

-.0021 

T= 

2.700 

E- 

-.0002 

E 1  = 

.0006 

£2  = 

-.0013 

T  = 

2.750 

E  = 

-.0002 

E  1  = 

.0005 

E2  = 

-.0016 

T  = 

2.800 

E  = 

-.0002 

£1  = 

.0005 

£2  = 

-.0014 

T  =  . 

2.850 

E  = 

-.0001 

£1  = 

.0004 

E2  = 

-.0012 

T= 

2.900 

b  = 

-.0001 

bl  = 

.0003 

E2  = 

-.0010 

T  = 

2.950 

£  = 

-.0001 

b  1  = 

.0003 

E2  = 

-.0009 

T  = 

3.000 

b  = 

-.0001 

bl  = 

.0003 

E2  = 

-.0007 

T = 

3.050 

E  = 

-.0001 

b  1  = 

.0002 

E2  = 

-.0006 

T  = 

3.100 

E- 

-.0001 

£1  = 

.0002 

£2  = 

-.0006 

T  = 

3.150 

E  = 

-.0001 

£  1  = 

.0002 

£2  = 

-.0005 

T  = 

3.200 

b  = 

-.0001 

tl  = 

.0001 

£2  = 

-.0004 

T = 

3.250 

.  E  = 

-.0000 

£1  = 

.0001 

£2  = 

-.0004 

T  = 

3.300 

£  = 

-.0000 

El  = 

.0001 

£2  = 

-.0003 
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Appendix  B 


ANALOG  COMPUTER  SOLUTION  OF  THIRD  ORDER  DIFFERENTIAL  EQUATION 
General  third  order  equation: 

E+aE+bE+cE  -  O 

a-  =  f,  +  '1  +  ri 

b  =  r,rz  +  r  r3  +  r2r3 

c  --  n  ri 


Scaling  factors: 


<xt  =  to 

cC  ^  —  oC  jr  *“  O .  / 


Basic  equation  for  an  integrator: 


Scaling  of  amplifier  stages: 

1 •  Amplifier  1 


•  • 


£ 


e+ 


+  C  t 


0.  = 


±n  QCt 

E  Q, 


«-*.  C+l 


t 


RzCf) 


C  -  aiaoct 

R»cft 


a 


IZ 


a 


is 


bPUO, 

oCt 

C  Cj, 

c<t 
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2 .  Amplifier  3 


E=/£dt 


3.  Amplifier  5 

E  =  /edt 

G 

&  =  =  '/ 

o  /  c 

I  =  a 


4 .  Amplifiers  21,  22,  23 

No  scaling  required  since  these  amplifiers  are  used  only  as  sign  changers. 


The  analog 


computer  schematic 


for  the  general  solution  is  shown  in  Fig. 


B-l. 
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Fig.  b'i  /Inaloj]  Computer  Set-up  for  Third  Order  Phase  Trajectories 


Appendix  C 


COMPUTER  SET-UP  AND  INITIAL  CONDITION  DERIVATION  FOR 
PHASE  TRAJECTORIES  IN  TRANSLATING  EIGENPLANE 

The  computer  set-up  as  used  on  the  Donner  3100  is  given  in  Fig.  C-l. 

The  initial  conditions  for  the  phase  trajectory  shown  in  Fig.  41 
were  for  a  point  in  the  complex  eigenplane  and  were: 

E  *  -  0.5,  E  »  +  1.0,  E  *+30 
o  o  o 

The  complex  eigenplane  was  translated  to  a  point  on  the  eigenvector 
corresponding  to : 

E  *  +4.0,  E  *  -  2.0,  E  *  1.0 

Therefore  the  initial  conditions  of  the  point  in  the  translated  eigen¬ 
plane  were: 


E  1 

*  E 

+  E 

-  +  3.5 

o 

o 

E  1 

*  E 

+  E 

*  -  1.0 

o 

o 

E  1 
0 

*  E 

+  E 

o 

=  +  4.0 

These  initial  conditions  gave  the  phase  trajectory  of  Fig.  42. 
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Fio.  C-l  /Analog  Computer  Set-up  Used  to  Obtain  Phase  Trajectory  of  Fu-».4Z 


Appendix  D 


DERIVATION  OF  THE  TRANCENDENTAL  EQUATION 
cos  3t  +  a  sin  j3t  38  e^a 


qe 


-(n+r^i 


+ 


rT  £ 


)(U-r>)  '  '  (r;-rO(r;-r5) 

This  expression  may  be  rewritten  as 


=  O 


A 


'  '\;t 


4 


Bf 


■(  r  ir,  ,•( 


where 


/l  = 


6  = 


'/L 


c  = 


-/? 


(rrr,)(rrr7)  (q-nX^-^)  <^-rzXr-r5) 


Dividing  through  by  C> 


/]_  -0^0*  5  <^>(  -fri+rs)t 

c  =« 


Here , 


_A_ _  r3(r-y;) 
c  >•;  o,-r> 


and 


S  _  n-  ^3 ) 


Let  r;  -  <ri_ 

rz  -  <*  +  j  (3 

B  -  <*  -  j  ^ 

Then  _  <aoC  -  o<X~  j  a  @  -  (=>z 

c  - 


O.OC  -  (3Z 

j  a  (Z 


wUr<L  CT  r:  *<*-<**- 3* 

^  (3 


(D-l) 


(D-2) 


(D-3) 


o- 
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Similarly,  p 


c r 


C  2  l 


J 


Multiplying  (D-3)  by  £ 


(^i  •  i)<- 


C 


C 


=  £ 


(D-4) 


Substituting  for  r^,  ,  r^,  ^  and  ^  in  (D-4), 


1  cr\.-J(3t  /|  <r  \  ie*  <«.-*.)  t 

~*/r  +(7+27)*  =  « 


z 


L(4>at+  £-*%*<-** 


2j 


cos  fit  -f  sin  fit  -  £ 
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thesD5674 

Er;ect0r  9e0metrv  Of  third  order  line 
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